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Similar  to  the  definition  of  tne  intensity  of  cavitation 
damage,  the  intensity  of  cavitation  bubble  collapse  is  defined 
as  tne  power  transmitted  per  unit  surface  area  of  the  babble 
when  the  collapse  pressure  is  a  maximum.  The  intensity  of 
collapse  is  given  by  tne  square  of  tr.e  maximum  collapse  pres¬ 
sure  divided  by  tne  acoustic  impedance  (p  C.  )  of  the  liquid,. 

■V  't 

Various  physical  effects  such  as  inertial  effects,  damping  due 
to  non-condensible  gases,  thermal  effects,  compressibility  ef¬ 
fects,  surface  tension  effects  and  viscous  effects  on  tne  maxi¬ 
mum  collapse  pressure  are  considered.  Tne  efficiency  of  damage 
given  by  the  Intensity  of  damage  divided  by  the  intensity  of 
collapse  is  snown  to  depend  principally  on  tne  dissolved  gas 
content  of  the  liquid  (using  tne  data  obtained  from  the  vibra¬ 
tory  cavitation  damage  apparatus).  It  is  also  snown  that  the 
vapor  itself  might  act  as  a  damper  near  boiling  point  since  the 
bubble  wall  temperatures  increase  rapidly  at  these  temperatures 
When  tne  vapor  pressure  becomes  important,  viscosity  also  seems 
to  affect  the  efficiency  of  damage  as  evidenced  oy  the  behavior 
of  aniline.  Tne  efficiency  of  damage  is  independent  of  the  sur 
face  tension  of  T'  liquid  within  tne  range  of  tests. 

Tne  above  results  from  tr.e  vibratory  experiments  are  used 
to  propose  a  modeling  technique  to  predict  the  r-  te  of  depth  of 
erosion  in  actual  operating  hydrodyr.ami c  systems.  One  of  the 
primary  questions  to  be  answered  is  the  magnitude  of  the  pres¬ 
sure  field  tha t  drives  the  bubble  to  collapse.  The  dependence 
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of  this  pressure  field  on  other  hydrodynamic  scaling  parameters 
such  as  Reynolds  number,  Mach  number  and  the  cavitation  parameter 
will  decide  the  success  of  the  proposed  model  technique. 

INTRODUCTION 

At  the  present  state  of  knowledge  it  is  possible  to  predict 
from  model  tests  the  various  hydrodynamic  forces  such  as  drag, 
lift,  thrust,  etc.  Similarly  it  would  be  highly  desirable  to 
predict  the  intensity  of . cavitation  damage  in  prototypes  by  con¬ 
ducting  model  tests  in  the  laboratory.  As  of  now  it  is  not  pos¬ 
sible  to  achieve  this  objective  because  there  are  no  scaling 
laws  that  relate  the  model-prototype  behavior.  It  is  the  pur¬ 
pose  of  this  report  to  formulate  such  scaling  laws  and  to  dis¬ 
cuss  the  physical  phenomena  that  are  scaled  by  these  laws. 

Figure  1  shows  the  flow  behind  a  circular  cylinder  and  the 
known  and  unknown  parameters  controlling  cavitation  damage.  The 
geometrical  and  kinematic  similarities  of  the  overall  cavity 
flow  are  controlled  by  the  cavitation  parameter  and  the  Reynolds 
number.  However  the  modeling  of  the  phenomenon  of  cavitation 
damage  requires  that  the  energy  of  collapse  of  individual  bub¬ 
bles,  the  transmission  of  the  energy  to  the  material  surface  and 
the  absorption  of  the  energy  by  the  material  in  its  deformation 
and  fracture  be  scaled  also.  Hence  the  problem  is  to  define  the 
above  phases  of  energy  transmission  and  absorption  quantitatively 
and  to  determine  the  parameters  that  control  the  efficiency  of 
this  process.  Figure  2.  The  intensity  of  material  damage  is 
defined  in  References  1  and  2  as 
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depth  of  erosion, 
time  of  erosion,  and 
erosion  strength. 


The  erosion  strength  is  defined  as  the  energy  absorbing  capacity 
of  the  material  per  unit  volume  under  the  action  of  the  erosive 
forces  (3).  If  a  similar  definition  for  the  intensity  of  bubble 
collapse  can  be  derived,  then  the  efficiency  of  this  process 
would  be  given  by 


I 

c 


[2] 


where  I  is  the  Intensity  of  bubble  collapse.  As  of  now,  such  a 
definition  for  the  intensity  of  bubble  collapse  does  not  exist. 
One  of  the  accomplishments  of  this  report  is  such  a  definition 
which  leads  to  logical  scaling  laws  for  the  proposed  technique 
of  modeling  cavitation  damage. 


INTENSITY  OF  BUBBLE  COLLAPSE 


It  is  assumed  in  this  report  that  the  collapse  pressure 
emanating  from  an  individual  transient  cavitation  bubble  causes 
the  erosion.  When  a  spherical  bubble  collapses  with  a  bubble 
wall  velocity  U  ,  the  collapse  pressure  as  derived  by  Rayleigh 
(4)  is: 
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is  the  maximum  collapse  pressure. 
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is  the  coefficient  of  compressibility  or  bulk 
modulus  of  liquid, 

is  the  density  of  liquid,  and 
is  the  bubble  wall  velocity. 
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wner-e  0.  is  the  sound  speed  in  liquid 
K/ 


The  compressibility  of  the  liquid  jb  would  be  greatly  re- 


■  ed  due  to  one  presence  of  cavitation  bubble; 


For  mis  reason 


the  sound  speed  of  the  liquid -vapor  bubble  mixture. 


wou  Id 


be  muen  less  than  0  .  The  actual  sound  speed  used  for  computing 

'G 

the  intensity  cf  collapse  should  be  C  ,  However,  very  little  is 

known  about  the  values  of  C  in  cavitation  clouds.,  Until  more 

m 

is  known  about  C  ,  the  sound  speed  in  liquid  Cl  will  be  used  for 

m  t 
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eomputat ions .  If  we  define  the  intensity  of  collapse  as  the 
radiated  power  per  unit  surface  area  of  the  bubble,  then 


since 


It  is  Interesting  to  note  that  the  Intensity  of  bubble  collapse 
is  similar  to  the  intensity  of  acoustic  power  radiated  from  a 
simple  source  (Figure  3).  The  dimensions  of  I  are  the  same  as 
I  which  are  in  power  per  unit  area. 

Q. 

The  above  discussion  shows  that  the  intensity  of  collapse 
varies  as  the  square  of  the  maximum  collapse  pressure  in  a 
given  liquid.  Hence,  it  is  important  to  consider  the  various 
physical  effects  that  control  this  maximum  collapse  pressure. 


Incorp orated 
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THE  MAXIMUM  COLLAPSE  PRESSURE 

The  maximum  collapse  pressure  depends  upon  the  shape  of 
the  bubble  and  the  distance  at  which  the  bubble  collapses. 
Furthermore,  the  growth  and  the  collapse  of  the  bubbles  depend 
upon  the  following  physical  phenomena  (5)>  (6): 

A.  Inertial  effects 

B.  Damping  effects 

C.  Thermal  effects 

D.  Compressibility  effects 

E.  Surface  tension  effects 

F.  Viscous  effects. 

Ir  order  to  quantitatively  evaluate  the  relative  influence  of 
tnese  effects  on  the  maximum  collapse  pressure,  we  will  assume 
that  the  Besant -Ray leigh  bubble  collapse  (5)  causes  the  damage. 
(It  is  recognized  that  the  Eisenberg  bubble  collapse  (7)>  (8), 
may  also  be  important.  Hence,  calculations  for  this  case  will 
be  useful. ) 

A,  Inertial  Effects 

Rayleigh  (4)  calculated  the  maximum  collapse  pressure  for 
a  spnerical  bubble  collapsing  in  an  infinite  liquid;  it  is 
given  by 


P 


[  53 
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where 

Pq  Is  the  pressure  at  infinity., 

R  is  the  maximum  bubble  radius  at  the 
start  of  collapse, 

and  R^  is  the  final  bubble  radius  at  the  end  of 

collapse  . 

This  maximum  collapse  pressure  occurs  at  a  distance  of  1.587  R^ 
from  the  center  of  collapse  (-4). 

B.  Damping  Effects 

Now  the  final  radius  to  which  the  bubble  collapses  depends 
upon  the  amount  of  non-condensible  gas  present  at  the  start  of 
the  collapse  and  on  whether  the  gas  is  compressed  isothermally 
or  adiabatically .  Assuming  isothermal  compression  of  the 
permanent  gas  inside  the  bubble,  Rayleigh  (4)  derives  the  fol¬ 
io-wing  relationship  between  the  final  collapse  radius  and  the 
initial  gas  content,  Qq: 
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Assumlng  adiabatic  compression,  Noltingk  and  Nepperas  (9)  de¬ 
rive  the  following  relationship: 
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where  k  is  the  ratio  of  specific  heats  for  the  permanent  gas  in¬ 
side  the  bubble.  The  value  of  k  for  air  is  A/3.  If  the  gas  is 
air,  then  for  adiabatic  compression. 
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Equations  [8]  and  [9]  are  graphically  shown  in  Figure  A.  It  is 
clear  that  the  maximum  collapse  pressure  and  hence  the  intensity 
of  collapse  depends  on  whether  the  bubble  collapse  is  isothermal 
or  adiabatic.  The  Plesset -Hsieh  criterion  (5)?  (10)  for  the 
collapse  of  bubbles  is  as  follows: 
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Isothermal 

Adiabatic 

where 

p  is  the  density, 

S  is  the  specific  heat  per  unit  mass  ,  and 
A  is  the  thermal  diffusion  length  . 

and  subscripts  l  and  g  stand  for  liquid  and  gas  respectively. 
Physically  speaking,  the  above  criteria  give  the  relationship  be¬ 
tween  the  heat  capacity  of  the  gas  inside  the  bubble  and  the  heat 
that  can  be  conducted  away  In  the  thermal  diffusion  layer  of  the 
liquid.  If  the  latter  is  very  large  compared  to  the  former,  the 
temperature  within  the  bubble  remains  constant.  If  it  is  very 
small  then  the  temperature  of  the  gas  inside  the  bubble  will  rise 
This  criterion  is  shown  in  Figure  5  for  various  test  liquids.  The 
thermodynamic  data  were  obtained  from  Reference  11  and  the  data 
for  argon  was  obtained  from  Reference  12.  At  20  kcs  frequency, 
bubbles  of  maximum  radii  of  a  millimeter  or  less  will  be  under 
isothermal  compression  while  collapsing,  for  all  the  liquids 
shown  in  Figure  5-  For  larger  bubble  sises,  the  collapse  would 
become  adiabatic.  For  this  reason,  the  model-prototype  correla¬ 
tion  should  take  into  account  the  Plesset -Hsieh  criteria.  How¬ 
ever  as  seen,  in  Figure  4,  adiabatic  collapse  produces  much  less 
pressures  and  may  not  be  important  as  far  as  cavitation  damage 
is  concerned. 
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Gas  content  of  the  bubble  -  The  partial  pressure  of  the 
gas  Inside  the  bubble  at  Its  maximum  radius.,  (in.  Equations  [8] 
and  [9])  consist  of  three  components : 

1,  The  gas  content  at  the  beginning  of  the  bubble 

growths  i.e.j  the  gas  inside  the  equilibrium 

nucleus  of  radius  R  . 

n 

2C  The  gas  dissolved  in  the  evaporated  liquid  during 
growth . 

3»  The  gas  that  diffuses  into  the  bubble  from  the 
surrounding  liquid  during  growth. 

Hence  the  partial  pressure  at  the  start  of  the  collapse,  Q  ,  will 
be  tne  sum  total  of  the  three  components  stated  above. 

Qq  =  Qi  +  Q2  +  Q3  [12] 


1.  Estimation  of  Qi  :  The  equation  of  static  equilib¬ 
rium.  for  spherical  bubbles  is  given  by 
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P, 


is  the  pressure  in  the  surrounding  liquid, 

is  the  equilibrium  radius  of  the  nucleus, 

is  the  surface  tension  of  the  liquid, 
is  the  vapor  pressure  of  the  liquid,  and 

Is  the  partial  pressure  of  the  gas. 
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Hence 


_  2.y 

P  +  ^  -  p 
o  R  v 
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Assuming  the  gas  inside  the  nucleus  to  be  a  perfect  gas. 


where 


p  %  ttR  3  =  Ni  0  T 
g  3  n  a 


Ni  is  the  number  of  moles  of  the  gas  , 

9  is  the  perfect  gas  law  constant  ,  and 

is  the  absolute  ambient  temperature  . 
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The  partial  pressure,  Qx,  exerted  by  these  molecules  when  the 
nucleus  grows  to  a  maximum  radius,  R  ,  is  given  by 


A  , 

Qx  ~  ttR  3  =  Ni  0  T 
3  m  a 


Qx  = 


P  +  — 1  -  p 

o  R  v 
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R 
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l  m 

C  13l 
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2.,  Estimation  of  Qa  :  The  mass  of  liquid  evaporated 
during  growth,  is  equal  to  the  mass  of  vapor  in  the  bubble  and 
is  g i v en  by 

4  „  3 

m,  =  r  f  R  p 

13  m 

where 

m  is  the  mass  of  liquid  evaporated 

/L/ 

is  the  density  of  vapor  inside  the  bubble 

The  mass  of  gas  which  was  originally  dissolved  in  this  mass  of 
liquid  before  evaporation  is 


m 

g 


a  m 
g  l 


4 

3 


7rR  3  p  a 
m  rv  g 


where  a  is  the  mass  fraction  of  gas  dissolved  in  the  liquid, 
g 

The  number  of  moles  of  gas  in  this  mass  is  given  by 


N2 
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M 


4  irR  3  p 
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M 

g 


where  M  is  the  molecular  weight  of  gas.  Similarly  the  number  of 
g 

moles  of  vapor  Nv  is  given  by 


N 
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4  3 

^ttR  p 

3  m  r 
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where  M.  is  the  molecular  weight  of  vapor. 
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Ns_ 

N 

v 


a  M 
g  v 


Assuming  the  perfect  gas  law  to  be  valid  both  for  gas  and  vapor 


where 


Qs  4  ttR  3  =  N3  0  T 

3  m  *  a 

4 

p  T  mR  3  =  N  0  T 
v  3  m  v  a 

Qs  is  the  partial  pressure  of  the  gas  exerted  by 
N3  gas  molecules 

p  is  the  vapor  pressure  inside  the  bubble 
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[14] 


3.  Estimation  of  Q3  :  While  the  bubble  grows  from  an 

equilibrium  nucleus  to  a  maximum  radius  R  ,  the  gas  dissolved  in 

m 

the  liquid  with  a  partial  pressure  of  Pq  will  diffuse  into  the 
bubble  from  a  gas  diffusion  layer  surrounding  the  bubble  and  is 
given  by  Bebchuck  (13)  as 
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wr.ere 


'D 
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( D  t)* 
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is  the  concentration  of  the  dissolved  gas, 
is  the  diffusivity  of  the  gas, 
is  the  diffusion  time,  and 
is  the  gas  diffusion  layer  thickness. 


Now  the  partial  pressure  of  the  gas  inside  the  bubble  at 
the  beginning  of  collapse  is  given  by  Equation  [  12]  which  becomes 
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Since  (R  /R  )  is  usually  small  compared  to  unity,  the  first  term, 
n*  m  ■ 

on  the  right  hand  side  of  the  above  equation  is  negligible  com¬ 
pared  to  the  other  two  terms.  Then 
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At  temperatures  much  lower  than  the  boiling  point  of  the  liquids 
the  vapor  pressure  is  small.  Hence 


o 


M 

v 

—  is  of  the  order  of  1 
g 

l. 

If  we  assume  that  the  gas  diffusion  layer,,  (D  t)2,  is  of  the 

g 

same  order  of  magnitude  as  R  then  the  Equation  [ l6]  may  be 
written  as  follows: 
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(D  t)  = 
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[  17l 


However  at  temperatures  close  to  the  boiling  point  of  the  liquids 
the  vapor  itself  might  act  as  a  damper  since  the  vapor  may  not 
fully  condense  during  collapse.  This  effect  will  be  discussed 
next  c 


C,  Thermal  Effects 

The  next  important  aspect  that  controls  the  collapse  of 
vapor  bubbles  is  the  heat  outflow  during  collapse.  The  heat  of 
condensation  is  given  in  (5)  as 
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4  3 

H  =  -  irR  3  p  L 
3  m  rv 


where 

Pv  is  the  density  of  vapor*  and 
L  is  the  latent  heat  of  vaporization. 

JL 

This  heat  flows  into  the  thermal  diffusion  layer  of  thickness  (Dt)2 
where  D  is  the  thermal  diffusivity  and  t  is  the  time  of  collapse. 
The  heat  balance  is  given  by 


R  L  p 

AT  =  — — — - 

3(Dt)%pt 

at  =  Rm  L  pv 

Tb  3  (Dtp  stPtrB 


[  18] 


where 

AT  Is  temperature  difference  between  the  liquid 
and  the  vapor* 

T  Is  the  boiling  point  of  the  liquid  at 
atmospheric  pressure*  and 

S  is  the  specific  heat  of  liquid  per  unit  mass. 

Assuming  a  millimeter  for  R^  and  20  p  seconds  for  growth  or 
collapse  time  in  various  liquids  tested  Figure  6  shows  that  the 
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relatlve  temperature  difference  increases  sharply  at  temperatures 
closer  to  the  boiling  points  of  the  liquids „  Thermodynamic  data 
are  obtained,  from  References  11  and  12. 

ti  that 

vi 

R  . 

the  bubble  wall,  i.e., 
the  liquid  adjacent  to 
the  vapor  adjacent  to 

Since  p_  ,/p  is  very  small,  it  is  evident  that  v  is  very  nearly 

V 

equal  to  P  for  most  practical  cases. 

Relat  ive  importance  of  dissolved  gas  cont  ent  and  vapor  pres  - 
sure  -  From  the  above  discussion,  it  Is  clear  that  when  AT/T 
Is  small,  the  damping  is  entirely  due  to  the  dissolved  gas  con¬ 
tent  In  the  liquid.  When  AT/T  is  increasing  with  ambient  temper- 

JD 

ature,  the  vapor  itself  may  act  as  a  damper.  At  these  tempera¬ 
tures  the  vapor  pressure  of  the  liquids  become  important. 

D,  Compressibility  Effects 

According  to  the  Rayleigh  bubble  collapse  mechanism,  the 
kinetic  energy  of  bubble  collapse  is  stored  in  the  bulk  compress¬ 
ibility  of  the  liquid  and  transmitted  back  to  the  material  in  the 


Furthermore  Plesset  (5)  has  shov< 
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where 


R  =  ~  Is  the  velocity  of 
the  interface 
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is  the  velocity  of 
the  interface,  and 

is  the  velocity  of 
the  interface. 
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form  of  a  short  range  shock.  In  order  to  Investigate  the  compre 
ibility  effects,  we  need  the  Mach  number  of  the  bubble  collapse. 
However,  there  are  three  sound  speeds  that  are  to  be  considered: 

1.  The  sound  speed  In  the  liquid 

2.  The  sound  speed  in  the  liquid -bubb le  mixture 

3 .  The  sound  speed  in  the  gas  within  the  bubble 

The  Mach  number  M  ,  may  be  defined  as 

w 


M 

w 


R 

C 


[191 


where  R  is  the  bubble  wall  velocity  and  C  is  any  of  the  sound 
speeds  that  is  Important.  If  we  again  consider  the  Rayleigh  bub 
ble  collapse,  then  R  is  given  by 


Since 


m 


»  i 


from  Equations  [6]  and  [7a] 
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For  isothermal  collapse 


[20] 


O  P  P 
_  _2  ° 
3  p  1h 


For  adiabatic  collapse 


During  the  analysis  of  experimental  results*  the  effect  of  these 
physical  parameters  will  have  to  be  considered. 

E.  Surface  Tension  Effects 

The  surface  tension  effects  may  be  scaled  by  the  Weber  num¬ 
ber*  W  *  for  the  motion  of  bubble  wall  given  by 
w 


p,  (R)3r„ 


w  = 

W 


[21] 


F.  Viscous  Effects 

Similarly  If  viscous  effects  were  to  become  important  in 
the  bubble  motions  during  growth  and  collapse*  then  the  corres¬ 
ponding  Reynolds  number,  R  *  for  the  bubble  motion  will  be  given 

w 


Rw  V 


[22] 


where  v  is  the  kinematic  viscosity  of  the  liquid. 
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EFFICIENCY  OF  CAVITATION  DAMAGE 


As  pointed  out  in  the  introduction,  the  efficiency  of  cavi¬ 
tation  damage  will  depend  upon  the  various  physical  effects  that 
control  the  maximum  collapse  pressure.  The  efficiency  of  cavi¬ 
tation  damage  may  be  defined  as 

i  S 


I 


h  = 


d 


P,C 


l  t 


From  the  previous  discussion, 
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[23] 


It  would  be  Interesting  to  see  which  non-dimensional  param¬ 
eters  in  Equation  [23]  really  influence  the  efficiency  of  cavi¬ 
tation  damage  for  a  typical  experimental  apparatus.  One  of  the 
most  widely  used  equipments  for  studies  on  cavitation  damage  is 
the  magnetostriction  vibratory  apparatus.  What  follows  is  a 
discussion  of  the  correlation  of  the  experimental  data  obtained 
with  this  apparatus. 
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CORRELATION  OF  EXPERIMENTAL  RESULTS  OF  VIBRATORY  TESTS 

The  basic  idea  in  these  experiments  is  to  vibrate  a  test 
specimen  in  a  given  liquid  contained  in  a  beaker  at  a  given 
amplitude  and  frequency  (see  Figure  7).  These  vibrations  are 
produced  in  most  cases  by  the  magnetostriction  oscillators. 

Since  1935  several  authors  (e.g.  References  14  through  19)  have 
conducted  these  tests  for  investigating  the  phenomenon  of  cavi¬ 
tation  damage.  The  basic  parameters  involved  in  these  tests  are 
the  beaker  dimensions.,  the  specimen  dimensions.,  amplitude  and 
frequency  of  vibrations,,  the  physical  and  chemical  properties  of 
the  material  and  liquid  used.  The  present  analysis  is  confined 
to  relatively  non-corrosive  liquids  and  metals.  Only  physical 
mechanisms  are  assumed  to  play  any  role. 

When  a  cylindrical  piston  vibrates  with  its  circular  face 
inside  the  liquid  as  shown  in  Figure  7a  the  maximum  pressure  Pq 
is  given  (Reference  20)  by 
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For  a  geometrically  similar  non-cavitating  system,  the  maximum 
pressure  depends  upon  the  frequency  parameter,  p,  given  by 


2ra 


where  a  is  the  radius  of  the  test  specimen  and  7.  is  the  wave 

AJ 

length  of  sound  in  the  liquid.  At  this  juncture,  it  Is  important 

* 

to  recognize  the  following  fact  pointed  out  by  Phillips.  It  is 
generally  observed  during  experiments  with  the  vibratory  appara¬ 
tus  that  there  is  a  rim  of  area  (at  the  periphery  of  the  eroded 
surface)  which  does  not  get  eroded  at  all  In  soft  aluminum  (see 
Figure  7b).  The  thickness  of  the  rim  is  indeed  of  the  order  of 
one-half  wave  length  as  pointed  out  by  Phillips  if  we  use  a  value 
of  about  100  fps  for  the  sound  speed  at  14000  cps  oscillations. 
This  magnitude  of  sound  speed  is  typical  for  air-water  mixtures. 
Here  again  the  importance  of  the  compressibility  of  the  two  phase 
mixtures  is  evident.  If  this  were  to  be  true,  then  the  measure¬ 
ment  of  rim  thickness  would  give  the  magnitude  of  sound  speeds  in 
cavitation  bubble  mixtures.  For  this  case  Equation  [23]  may  be 
written  as 
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[25] 


-* 

O.M.  Phillips,  personal  discussion  with  the  author. 
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where 


and 


r 
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is  the  volume  concentration  of  dissolved 
gas  (damping  due  to  gas  content) 
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motion 

relative  vapor  pressure 


relative  diffusion  length 

relative  -bubble -wall- temperature  -in- 
during  collapse  (thermal  effects) 

Mach  number  of  the  bubble  wall  motion 
ssibility  effects) 

Weber  number  of  the  bubble  wall 
(surface  tension  effects) 

Reynolds  number  of  the  bubble  wall 
(viscous  effects) 


is  the  frequency  parameter. 


The  experimental  results  of  several  authors  are  shown  in 
Figure  8  along  with  the  range  of  variables  involved  in  these  tests. 
The  rate  of  weight  loss  is  plotted  against  temperature  in  the 
melting -boiling  range  for  each  liquid.  All  these  experiments  were 
conducted  at  atmospheric  pressure  except  for  one  case  in  which 
Kerr  and  Leith  (l6)  conducted  experiments  at  2.4  atmospheres  ab¬ 
solute.  The  table  attached  to  this  figure  gives  all  the  avail¬ 
able  information  on  the  liquid,,  material,  amplitude,  frequency 
and  pressure  used  in  these  tests.  These  data  provide  substantial 
information  to  understand  the  relative  importance  of  the  various 
parameters  contained  in  Equation  [25].  The  surface  tensions  of 
these  liquids  vary  from  25  dynes/cm  (for  aniline)  to  75  dynes/cm 
(for  water),  the  viscosity  from  0.3  centipoises  to  10  centipoises. 
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the  sound  speed  from  1320  meters /sec  to  1650  meters/sec  and  the 
density  from.  0.8  gm/cm3  to  1.0  gm/cm'3  .  These  data  were  obtained 
from  the  International  Critical  Tables  and  from,  the  Handbook  of 
Physics  and  Chemistry. 

The  efficiency  as  given  in.  Equation  [25]  was  calculated 
from  these  data  and  plotted  against  the  dissolved  air  content  as 
shown  in  Figure  9.  The  rate  of  depth  of  erosion  was  obtained 
from,  rate  of  weight  loss  divided  by  the  density  of  the  material 
and  the  area  of  specimen.  The  rate  of  erosion  is  time  dependent 
(l).,  (21),  (22),  Our  experimental  data  were  obtained  in  the 
steady  state  whereas  all  the  other  investigators  did  not  separate 
the  time  effect.  This  Is  an  obvious  limitation.  The  erosion 
strength  values  for  the  metals  have  been  obtained  from,  previously 
published  results.  In  some  cases  as  indicated  in  Figure  9,  an 
equivalent  erosion  strength  has  been  used  to  take  into  account 
tne  time  effect.  For  water  the  dissolved  air  content  was  measured 
with,  the  Van  Slyke  apparatus.  The  dissolved  air  content  decreases 
with  testing  time  at  a  constant  temperature  ir:  the  vibratory  rests 
as  shown  in  Figure  10  and  reaches  a  steady  value.  The  steady 
value  depends  upon  the  test  temperature  as  shown  in  Figure  11  „ 

In  Figure  10,  the  published  values  of  air-  content  in  water1  at 
various  temperatures  are  also  shown  for  comparison.  The  air  con¬ 
tent  for  aniline,  toluene  and  benzene  are  given  in  International 
Critical  Tables  for  room  temperature,,  The  air  content  was  as  - 
sumed  to  vary  inversely  as  the  temperature  and  was  calculated  for 
every  other  temperature  as  shewn,  in  Figure  11. 
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The  rate  of  damage  Increases  with  temperature  and  then  de¬ 
creases  with  increasing  temperature  as  shown  in  Figure  8=  The 
peak  damage  occurs  around  50  to  70  percent  of  the  boiling  point. 
As  shown  in  Figure  6,  in  this  temperature  range  the  relative 

bubble  wall  temperature,  AT/T_,  seems  to  increase  rapidly  with 

B 

the  result  that  the  vapor  pressure  also  becomes  important.  All 
the  data  shown  In  Figure  9  are  for  the  case  when  this  effect  is 
not  Important  (i.e„,  when  the  ambient  temperature  is  below  the 
temperature  that  produces  the  peak  damage) „  Above  this  tempera¬ 
ture,  the  relative  vapor  pressure  becomes  the  important  corre¬ 
lating  parameter  as  shown  in  Figure  12  in  which  the  values  of 

r\/r\  are  plotted  against  p  /P  for  water,  benzene,  toluene 
peak  v  o 

and  aniline.  Aniline  seems  to  behave  differently  as  compared  to 
the  other  three  liquids.  One  possible  explanation  is  that  the 
viscosity  of  aniline  is  one  order  of  magnitude  higher  than  that 
for  the  other  three  liquids.  Viscous  damping  seems  to  be  Impor¬ 
tant  at  higher  temperatures  when  vapor  pressure  plays  a  dominant 
role.  However  at  lower  temperatures  the  viscosity  seems  to  be 
unimportant  as  shown  by  the  correlation  in  Figure  9.  This  is 
further  supported  by  the  experiments  with  solutions  of  non- 
Newtonian  additives  shown  in  Figure  13.  Vibratory  cavitation 
damage  tests  were  conducted  in  solutions  of  sodium  carboxy- 
methy lcellulose .  The  concentrations  of  these  additives  were 
varied  such  that  the  viscosity  of  these  solutions  would  vary  ten¬ 
fold.  No  noticeable  change  in  rate  of  damage  were  observed  at 
room  temperature. 
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Considering  the  fact  that  the  results  presented  in  figure  9 
belong  to  a  wide  variety  of  data  colie. tea  by  different  authors 
using  different  experimental  conditions,  the  efficiency  of  damage 
seems  to  depend  primarily  on  dissolved  gas  content  of  the  liquid 
at  temperatures  well  below  the  boiling  point.  Additional  infor¬ 
mation  presented  in  Figure  12  seems  to  Indicate  that  the  relative 
vapor  pressure  Is  also  important  at  temperatures  close  to  boiling 
point  The  characteristic  behavior  of  aniline  indicates  the 
importance  of  viscosity  at  higher  temperatures.  Additional  ques¬ 
tions  as  to  whether  the  compressibility  and  surface  tension  ef¬ 
fects  are  Important  or  not  remains  to  be  verified  by  additional 
experiments  since  the  range  of  values  of  surface  tension  and 
sound  speeds  in  the  present  case  considered  are  not  wide  enough . 
Fortner more,  the  gas  content  in  the  liquids  must  be  actually 
measured  while  other  properties  are  accurately  controlled.  Hence 
further  investigations  to  refine  these  correlations  and  to  check 
trie  significance  of  compressibility  and  surface  tension  effect? 
are  necessary. 

However  the  correlations  shown  in  Figures  9  and  12  already 
explain,  a  priori,  a  few  of  the  significant  experimental  results 
that  neeaed  clarification; 

1,  For  example,  the  variation  of  damage  intensity 
with  the  square  of  the  amplitude  of  motion  ((21)  and  (23))  can. 
be  explained  as  follows.  The  intensity  should  be  proportional 
to  the  square  of  the  amplitude  according  to  Equation  [  25]  if  ail 
tree  terms  on  the  right  hand  side  of  this  equation  remain  constant 
ar.d  independent  of  the  amplitude. 
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2.  In  hydrocarbon  liquids  such  as  benzene,  toluene, 
aniline,  etc.,  the  intensity  of  damage  Is  generally  observed  to 
be  much  lower  as  compared  to  water.  This  can  be  explained  in 
terms  of  the  higher  solubility  of  air  in  hydrocarbon  liquids. 

3.  Wilson  and  Graham  (18)  observed  that  the  rate  of 
damage  strongly  depended  on  the  value  of  p, C  (acoustic  impedance 
of  the  liquid)  and  an  explanation  of  this  effect  is  obvious  from 
the  present  analysis.  The  non-linearity  in  their  correlation  may 
be  explained  in  terms  of  the  variation  in  solubility  of  air  in  the 
various  liquids  they  studied. 

7.  Another  interesting  observation  by  Wilson  and 
Graham  (l8)  is  that  the  rate  of  damage  did  not  vary  with  viscosity 
over  a  wide  range  in  their  experiments  with  water  and  glycerene 
mixtures  at  room  temperature.  This  confirms  our  own  experiments 
with  non -Newtonian  additives  shown  in  Figure  13.  These  results 
verify  that  viscous  effects  are  not  important  at  lower  tempera¬ 
tures  . 

5.  The  dependence  of  damage  on  temperature  is  gener¬ 
ally  explained  (see  for  example  (17)  and  (19))  as  follows.  At 
lower  temperatures  the  damage  increases  with  temperature  because 
the  solubility  of  air  decreases  with  temperature.  However  at 
higher  temperatures  the  vapor  pressure  becomes  Important  in 
dampening  the  collapse  of  bubbles.  These  explanations  are  quanti¬ 
tatively  demonstrated  in  Figures  9  and  12. 
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If  we  neglect  the  compressibility  aria  surface  tension  ef¬ 
fects ,  we  may  rewrite  Equal  lor  [23]  as 


1  D 


P„ 
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=  F  T  ,  ~  ,  ft,  ,  Flow  Parameters 

P  2  g  Po  w 
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1  "i 


When  vapor  pressure  is  rot  important. 


p  = 


F  ,  Flow  Parameter, 


u  27] 


HYDRODYNAMIC  CAVITATION  DAMAGE 


Now  tne  question  Is  how  to  maKe  use  of  the  correlations  ob¬ 
tained  for  the  vibratory  tests  to  derive  similar  relationships 
governing  the  efficiency  of  cavitation  damage  in  actual  hydro¬ 
dynamic  flew  systems..  Equation  [27]  may  be  rewritten  for  this 
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'is  the  cavitation  parameter  ar..d 
is  tne  Reynolds  number  of  tne  flow „ 
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How  ever  the  magnitude  of  the  characteristic  pressure  field  Pq 
that  forces  the  bubble  to  collapse  is  not  readily  known.  Let  us 
consider  a  practical  example  shown  in  Figure  1  where  cavitation 
damage  is  produced  in  the  wake  of  a  circular  cylinder.  The 
Reynolds  number  and  cavitation  parameter  will  be  given  by 

V  D 
00 

R  =  - 

e  v 


and 


a  = 


The  turbulent  pressure  field  would  be  given  by 


P 

turbulent 


u 


[29] 


Similarly  the  stagnation  pressure  also  will  be  given  by 


P 

o 


stagnation 


[30] 


where 

V  is  the  free  stream  velocity 
00  0 

p^  is  the  free  stream  pressure 

u1  is  the  mean  turbulent  velocity  fluctuation,  and 
D  is  the  diameter  of  the  cylinder. 
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If  either  of  the  above  two  relationships  were  to  be  true.,  the 
intensity  of  damage  should  vary  as  the  fourth  power  of  the  free 
stream  velocity  (see  Figure  14),  However,  Knapp  (24)  found  that 
the  intensity  of  damage  varied  as  the  6th  power  of  the  free 
stream  velocity.  The  value  of  the  exponent  is  known  to  vary 
from.  5  to  7  according  to  Kerr  and  Rosenberg  (25),  Thiruvengadatr 
(26)  found  that  the  exponent  varied  from  4,5  to  8.5  at  a  const an 
value  of  the  cavitation  parameter  depending  upon  the  eavitatlng 


These  considerations  lead  us  to  conclude  that  an  entirely 
different  approach  may  be  needed  to  explain  the  high  power 
dependency  on  velocity.  Lightniil  (27,  28)  analyzed  the  mech¬ 
anism  of  conversion  of  energy  from  kinetic  energy  of  fluctuating 
shearing  motions  (Reynolds  stresses)  into  the  acoustic  energy  of 
fluctuating  longitudinal  motions  and  derived  the  following  equa¬ 
tion  from  dimensional  analysis: 


I 

a 


II 

Ct 


[31] 


where  I  is  the  acoustic  intensity  radiated  by  the  turbulent 

Q. 

eddies  acting  as  quadrupole  sources.  If  we  assume  that  the 
cavitation  bubbles  get  entrained  by  these  eddies,  then  the  eddie 
would  correspond  to  oscillators  driving  the  bubbles  to  collapse 
(Figure  14).  Then  just  as  in  vibratory  tests  the  efficiency  of 


damage  may  be  written  as 
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where  M  Is  the  Mach  number  of  the  flow,  There  are  no  systematic 
experimental  data  available  to  check  the  validity  of  the  above 
relationships . 


Furthermore  cavitation  damage  also  depends  strongly  on  the 
value  of  the  cavitation  parameter  at  a  given  velocity  (29,  26). 
Intensity  of  damage  is  maximum  f or  a  specific  a  value  at  any 
given,  velocity.  These  considerations  indicate  that  two  sets  of 
curves  are  necessary  in  order  to  predict  a  model  prototype  rela¬ 
tionship  as  shown  in  Figure  15 .  The  proposed  method  of  model 
testing  needs  verification. 


CONCLUSIONS 


The  following  conclusions  may  be  drawn,  from,  the  previous 
analyses : 

1  The  Intensity  of  bubble  collapse  is  defined  as  the 
power  radiated  per  unit  surface  area  of  the  bubble,.  It  Is  given 
by  the  square  of  the  collapse  pressure  divided  by  the  acoustic 
impedance  of  the  liquid. 

2.  The  efficiency  of  damage  Is  defined  as  the  ratio 
of  the  intensity  of  cavitation  damage  to  the  intensity  of  bubble 
collapse . 
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3 .  In  the  vibratory  tests*  the  efficiency  of  damage 
is  primarily  controlled  by  the  isothermal  damping  of  the  dis¬ 
solved  non-condensible  gases  at  lower  ambient  temperature  of 
the  liquid.  At  higher  temperatures  closer  to  the  boiling  point* 
the  vapor  of  the  liquid  itself  plays  an  important  role.  This  is 
because  the  relative  bubble  wall  temperatures  increase  rapidly 
at  these  temperatures. 

4.  When  the  vapor  pressure  is  important*  the  viscosity 
of  the  fluid  also  seems  to  be  important  as  shown  by  the  behavior 
of  aniline, 

5.  Within  the  range  of  tests*  the  efficiency  of  dam¬ 
age  seems  to  be  independent  of  the  surface  tension  of  the  liquid. 

6.  The  above  results  of  the  vibratory  tests  have  been 
extrapolated  to  propose  a  modeling  technique  to  predict  the  in¬ 
tensity  of  cavitation  damage  in  actual  hydrodynamic  systems. 

The  proposed  modeling  technique  needs  experimental  verification. 
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FIGURE  2  -  PHASES  OF  ENERGY  TRANSMISSION 


P2  p  .  -  LIQUID  DENSITY 

INTENSITY  OF  COLLAPSE  I  =  — -  £ 

c  p£Cf  C  -  SOUND  SPEED  IN  LIQUID 


FIGURE  3  -  DEFINITION  OF  INTENSITY  OF  BUBBLE  COLLAPSE 
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FIGURE  4  -  BUBBLE  COLLAPSE  PRESSURES  AS  A  FUNCTION  OF  GAS  CONTENT  BOTH  FOR 
ISOTHERMAL  AND  FOR  ADIABATIC  DAMPING 


FIGURE  5  -  PLESSET-HSIEH  CRITERION  FOR  VARIOUS  LIQUIDS  TESTED 
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FIGURE  6  -  THERMAL  EFFECTS  ON  VAPOR  BUBBLE  GROWTH  AND  COLLAPSE 
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FIGURE  7-  PARAMETERS  INVOLVED  IN  MODELLING  VIBRATORY  CAVITATION 
DAMAGE  EXPERIMENTS 


HYDRO  NAUTICS,  INCORPORATED 


pro 

—jfcpJiA  -  \ - -I - \ 

m  /'N^Ut  ^ 


6/ - -pi 


AMBIENT  TEMPERATURE  -  °C 


INVESTIGATOR 


TE  ST  IFREOUENCY  I  AMPLITUDE 


SCHUM6,  PETERS  ALUMINUM 
[AND  MILLIGAM  51  -  ST  1  WA™ 


MATERIAL  I  LIQUID 


DATA  OBTAINED  BY  MRS.  5.  W  WHITE  WITH  THE 
HYDRONAL TICS'  MAGNETOSTRICTION  APPARATUS 


LIQUID  PROPERTY  RANGE  COVERED  BY  CAVITATION  DAMAGE 


VISCOSITY  |  SUR-ACE  DENSITY 


VISCOSITY 

Cenfi  poises 


SURFACE 
TENSION 
dynes  'cm 


WATER 


ANILINE 

BENZENE 

TOLUENE 


FIGURE  8  -  SUMMARY  OF  RESULTS  OF  VIBRATORY  TESTS 
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FIGURE  9  -  CORRELATION  OF  EFFICIENCY  OF  DAMAGE  WITH  AIR  CONCENTRATION 
IN  VARIOUS  LIQUIDS 
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FIGURE  10  -  VARIATION  OF  DISSOLVED  AIR  CONTENT  WITH  TEST  DURATION  IN 
VIBRATORY  TEST 


VOLUME  OF  AIR  ATS. 
VOLUME  OF  SOLVENT 


TOLUENE 


BENZENE 


O  WATER  (37th  Ed.  Chem.  and  Phys.  HDBK  Values) 
0  WATER  (Measured  Values  after  Cavitation) 


AMBIENT  TEMPERATURE  -  °C 
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FIGURE  12  -  CORRELATION  OF  EFFICIENCY  WITH  VAPOR  PRESSURE  AT  TEMPERATURES 
CLOSE  TO  BOILING  POINT _ — - 
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FIGURE  13  -  EFFECT  OF  NON-NEWTONIAN  ADDITIVE  ON  CAVITATION  DAMAGE  RATE 
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FIGURE  14  -  HYDRODYNAMIC  CAVITATION  DAMAGE 
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FIGURE  15  -  PROPOSED  METHOD  OF  MODELING 
CAVITATION  DAMAGE 
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